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C/^ . ABSTRACT 

Aims. The EROS -2 project has been designed to search for microlensing events towards any dense stellar field. The densest parts of 

I I the Galactic spiral arms have been monitored to maximize the microlensing signal expected from the stars of the Galactic disk and 

^ , bulge. 

' Methods. 12.9 million stars have been monitored during 7 seasons towards 4 directions in the Galactic plane, away from 

\ the Galactic center. 

i ■ Results. A total of 27 microlensing event candidates have been found. Estimates of the optical depths from the 22 best 

' events are provided. A first order interpretation shows that simple Galactic models with a standard disk and an elongated bulge are 

I in agreement with our observations. We find that the average microlensing optical depth towards the complete EROS -cataloged stars 

^\ , of the spiral arms is t = 0.51 + .13 x 10~^, a number that is stable when the selection criteria are moderately varied. As the EROS 

■ catalog is almost complete up to Ic = 18.5, the optical depth estimated for the sub-sample of bright target stars with Ic < 18.5 

^ ■ (t = 0.39 + . 1 1 X 10"^) is easier to interpret. 

, Conclusions. The set of microlensing events that we have observed is consistent with a simple Galactic model. A more pre- 

^ j ■ cise interpretation would require either a better knowledge of the distance distribution of the target stars, or a simulation based on a 

' Galactic model. For this purpose, we define and discuss the concept of optical depth for a given catalog or for a limiting magnitude. 
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1. Introduction 

After the first reports of microlensing candidates 
(Aubourg etal 1993, Alcock etal 1993, Udsihki etal 1993), 
the EROS team has performed extensive microlensing surveys 
from 1996 to 2003, that monitored the Magellanic clouds and 
large regions in the Galactic plane. The EROS-2 search for 
lensing towards the Magellanic clouds (Tisserand et ah 2007) 



" Now at LPNHE, 4 place Jussieu, F-75252 Paris Cedex 5, France 
** Now at San Pedro de Atacama Celestial Exploration, Casilla 21, 
San Pedro de Atacama, Chile 
Deceased 

Now at Mount Stromlo Observatory, Weston P.O., ACT, 2611, 
Australia 



2 



EROS collaboration: Microlensing towards the spiral arms 



yielded significant upper limits on the fraction of the Milky Way 
halo that can be comprised of dark objects with masses between 
10"^ Mo and IOMq. For objects of mass OAMq the 95% CL 
limit is 8%, in conflict with the suggestion by the MACHO 
coUaboration (Alcock et al 2000) that between 7% and 50% 
of the halo is made up of such objects. The EROS-2 search for 
microlensing of Galactic Bulge clump giants yielded 120 events 
(Hamadache et al 2006) giving a Galactic-latitude dependent 
optical depth of 



r/lO"^ = (1.62 ± 0.23)exp[-a(|Z?| - 3°)] , 
with 

a = (0.43 ±0.16) deg-^ . 



(1) 



(2) 



This optical depth agrees with Galactic models 
(Evans & Belokurov 2002 ; Bissantz et al 1997) and with 
the resuhs of the MACHO (Popowski et al 2005) and Ogle-II 
(Sumi et al. 2006) collaborations. The duration distribution of 
the events discovered by the three collaborations have been 
recently analyzed by Calchi Novati et al 2008 to constrain the 
Galactic Bulge Initial Mass Function. 

Our team has devoted about 15% of the observing time dur- 
ing 7 seasons to the search for microlensing events towards 
the Galactic Spiral Arms (GSA), as far as 55 degrees in lon- 
gitude away from the Galactic center. In our previous publica- 
tions (Derue et al. 1999, Derue et al. 2001, hereafter referred as 
papers I and II) describing the detection of respectively 3 and 7 
events, our attention was called on a possible optical depth asym- 
metry, accompanied by an asymmetric event dynamics with re- 
spect to the Galactic center. This marginal efl'ect (a 9% proba- 
bility to be accidental) could be interpreted as an indication of 
a long Galactic bar within the bulge. Its investigation required a 
significant increase in the number of events. 

In addition to the observing time increase (by more than a 
factor 2), we improved our catalog of monitored stars by increas- 
ing the limiting magnitude as well as by recovering some fields 
and sub-fields that were not analyzed previously. These improve- 
ments allowed us to recover another factor ~ 1.5 in sensitiv- 
ity. Moreover the discrimination power for microlensing event 
identification has been significantly improved, partly because the 
light curves are longer and thus provide a better rejection of re- 
current variable objects. 

A specific difficulty in the analysis of the spiral arms survey 
comes from the poor knownledge of the source distance distri- 
bution; in contrast with the LMC, the SMC and the Galactic cen- 
ter red giant clump, the monitored sources in the Galactic disk 
span a wide range of distances (±5 kpc according to preliminary 
studies, see Sect. 6.4.2). Their mean distance is also uncertain 
and has been estimated tohcl ±lkpc (Derue 1999b). We define 
in this paper the notion of "catalog optical depth" (Sect. 9) and 
provide all the necessary data to test Galactic models. 

2. Microlensing basics 

Gravitational microlensing (Paczynski 1986) occurs when a 
massive compact object passes close enough to the line of sight 
of a star, temporarily magnifying the received light. In the ap- 
proximation of a single point-like object acting as a deflector on 
a single point-like source, the total magnification of the source 
luminosity at a given time t is the sum of the contributions of 
two images, given by 



A(t) = 



u{tf + 2 

u{t) yju(t)^ + 4 ' 



(3) 



where u(t) is the distance of the deflecting object to the unde- 
flected line of sight, expressed in units of the "Einstein Radius" 
Re: 



Re 



-4 



4GM 



Lx(l - x) , 



(4) 



4.54 A.U. X 



M 



10 kpc 



[x(l-x)]-2 
0.5 



Here G is the Newtonian gravitational constant, L is the distance 
of the observer to the source and xL is its distance to the de- 
flector of mass M. The motion of the deflector relative to the 
line of sight makes the magnification vary with time. Assuming 
a deflector moving at a constant relative transverse speed Vt, 
reaching its minimum distance uq (impact parameter) to the un- 
deflected line of sight at time to, u{t) is given by 



u{t) = ^|ul + {{t-tQ)|tEf, 



(5) 



where tE = the "lensing time scale", is the only measurable 
parameter bringing useful information on the lens configuration 
in the approximation of simple microlensing: 



tE{days) = 79. 



Vt 



100 km/ s 



M 



Mo 



lOkpc 



[x(l-x)]-2 
0.5 



.(6) 



This simple microlensing description can be broken in 
many diff'erent ways : double lens (Mao & Stefano 1995), ex- 
tended source, deviations from a uniform motion due ei- 
ther to the rotation of the Earth around the Sun (parallax 
effect)(Gould 1992, Hardy & Walker 1995), or to the orbital mo- 
tion of the source around the center-of-mass of a multiple sys- 
tem, or to a similar motion of the deflector (see for example 
MoUerach & Roulet 2002). 

The optical depth r towards a particular set of target stars is 
defined as the average probability for the line of sight to inter- 
cept the Einstein disk of a deflector (magnification A > 1.34). 
This probability is independent of the deflector mass function, 
since the surface of the Einstein disk is proportional to the de- 
flector's mass. When the target consists of a population of stars, 
the measured optical depth is obtained from 



1 



T = 



TT ^ tE 

"9 Zj 7(t^' 



No,s^To,s2^^^^e{tE) 



(7) 



where A^^^^ is the number of monitored stars; IsTobs is the du- 
ration of the observing period; e{tE) is the average detection ef- 
ficiency of microlensing events with a time scale tE, defined as 
the ratio of detected events to the number of events with < I 
whose magnification reaches its maximum during the observing 
period. Similarly, the event rate corrected for the detection effi- 
ciency is 



1 



Nobs^Tobs e{tE) 



(8) 



3. Experimental setup and observations 

The telescope, the camera and the observations, as well as the 
operations and data reduction are described in paper I and ref- 
erences therein. The star population locations and the amount of 
data collected towards the 29 fields that have been monitored in 
four diff'erent regions (j3 Set, y Set, y Nor and 6 Mus) are given in 
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Fig. 1 and table 1. Taking into account the dead zones, the lower 
efficiency sectors of our CCDs and the blind zones around the 
brightest stars, we estimate that 75 ± 4% of the total CCD area 
(0.95 deg^) was effectively sensitive. This number was obtained 
by estimating the excess of 10 x 10 pixel domains (6'' x 6'') 
containing zero star, with respect to the number of void domains 
expected from the Poissonian distribution of the stellar number 
density. It is in agreement with the ratio between the total num- 
ber of detected stars (summed over all fields) and the number 
extrapolated from the stellar density observed in the CCD best 
zones. We took exposures of 120 ^ towards JS Set, y Set and y Nor 
and 180 5 towards Mus. The observations span a period of 
ATobs = 2325 days, starting July 1996 and ending October 2002; 
369 measurements per field were obtained on average in each 
of the Reros and Beros bands. Our fields were calibrated us- 
ing the DENIS catalog (Epchtein et al 1999) and the calibration 
was checked with the OGLE-II catalog (Udalski et al 2000b). 
We found that Reros and Beros bands are related to the Cousins 
I and Johnson V magnitudes through the following color equa- 
tions, to a precision of ~ 0.1 mag: 



R 



EROS 



Ic, Beros = Vj - OA(Vj - Ic). 



(9) 



Figure 2 shows the observation time span and the average sam- 
pling for the four diff'erent directions. 



Table 1. Characteristics of the 29 fields which were monitored in the 
EROS spiral arm program: Locations of the field centers, average sam- 
pling (number of photometric measurements per light curve and per 
color) and number of stars monitored for each field. The observing time 
was ATobs = 2325 days. The total numbers of observed stars towards 
7 Nor and 6 Mus are smaller than the sum of the numbers given for 
each field because of some overlap between contiguous fields. The total 
fields of view (f.o.v.) are the areas effectively monitored (0.71 deg^ per 
field, see text), corrected for the overlap between fields. 



Field 


a° (/2000) 


S° (/2000) 




/° 


Nmeas 


Notsim 




Set Exposure 


= 120s. f.o.v. 


= 4.3deg^. 




268 


3.00 


bs300 


280.8417 


-7.6814 


-1.75 


25.20 


269 


0.48 


bs301 


280.8625 


-6.2283 


-1.11 


26.51 


266 


0.47 


bs302 


281.5667 


-7.3792 


-2.25 


26.80 


272 


0.50 


bs303 


281.5833 


-5.9264 


-1.60 


27.09 


261 


0.47 


bs304 


282.3375 


-6.7642 


-2.70 


26.71 


269 


0.52 


bs305 


283.1083 


-6.5956 


-3.26 


27.19 


271 


0.54 


7 Set Exposure 


= 120s. f.o.v. 


= 3.6 deg^ 




277 


2.38 


gs200 


277.0125 


-14.8517 


-1.64 


17.72 


282 


0.47 


gs201 


277.8125 


-14.2439 


-2.12 


18.00 


266 


0.49 


gs202 


277.8875 


-12.8147 


-1.52 


19.30 


291 


0.49 


gs203 


278.5917 


-14.5275 


-2.92 


18.09 


281 


0.46 


gs204 


278.6167 


-13.0753 


-2.28 


19.40 


267 


0.47 


y Nor Exposure 


= 120s. f.o.v. 


= 8.4 deg^ 


454 


5.24 


gn400 


242.4375 


-53.1175 


-1.17 


330.49 


496 


0.42 


gn401 


244.5917 


-51.7453 


-0.99 


332.04 


475 


0.41 


gn402 


243.7375 


-53.0764 


-1.59 


330.74 


463 


0.45 


gn403 


245.6167 


-52.1056 


-1.69 


332.24 


420 


0.42 


gn404 


244.7875 


-53.4439 


-2.29 


330.94 


435 


0.43 


gn405 


246.7167 


-52.3506 


-2.35 


332.54 


445 


0.44 


gn406 


245.9750 


-53.7314 


-2.99 


331.23 


443 


0.46 


gn407 


247.8792 


-52.4789 


-2.95 


332.93 


443 


0.47 


gn408 


247.1750 


-53.8661 


-3.60 


331.63 


453 


0.47 


gn409 


243.9625 


-54.8125 


-2.86 


329.82 


482 


0.47 


gn410 


245.1250 


-55.0717 


-3.59 


329.93 


443 


0.46 


gn411 


242.4042 


-55.1686 


-2.54 


328.78 


449 


0.48 


e 


Mus Exposure 


= 180s. f.o.v. 


= 3.8 deg^ 




375 


2.28 


tmSOO 


201.7667 


-63.0383 


-0.47 


306.98 


391 


0.44 


tmSOl 


202.8250 


-63.5781 


-1.07 


307.37 


355 


0.44 


tm502 


203.7167 


-64.1750 


-1.72 


307.66 


376 


0.47 


tm503 


200.9917 


-64.9978 


-2.36 


306.38 


375 


0.36 


tm504 


198.0500 


-64.1136 


-1.35 


305.22 


392 


0.43 


tm505 


199.0625 


-64.6806 


-1.96 


305.60 


360 


0.43 


Total 


12.9 
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Fig. 2. Time sampling for each monitored direction: weekly average 
number of measurements per star since January Irst, 1996. 



4. The catalogs 

The catalogs of monitored stars have been produced following 
the procedure described in papers I and II, based on the PEIDA 
photometric software (Ansari 1996). All objects are well identi- 
fied in both colors and unambiguously associated between these 
two colors. We have removed objects that suff'er from a strong 
contamination by a nearby bright star; the contribution to the 
background flux from such a nearby star at the position of the 
object should not exceed 150% of its peak flux. 

The seven season data set contains 12.9 million objects mea- 
sured in the two colors: 3.0 towards J3 Set, 2.4 towards y Set, 
5.2 towards y Nor and 2.3 towards 6 Mus. The number of moni- 
tored stars was increased by ~ 50% since the analysis of papers 
I and II, by producing a richer catalog from a wider choice of 
good quality images than available before. We were also able to 
solve some technical problems that prevented us from producing 
the catalog for some fields (Tisserand 2004, Rahal 2003). The 
recovered stars are mainly faint stars with a comparatively low 
microlensing sensitivity. 



4.1. Completeness, blending 

We have compared a subset of the gs201 EROS field cata- 
log (Fig. 3a) with the catalog extracted from the deeper HST- 
WFPC2 (Wide Field Planetary Camera 2) images (Fig. 3b) 
named U6FQ1102B (exposure 210s with filter F606W) and 
U6FQ1104B (exposure 126s with filter F814W), centered at 
(a = 277.6281°, S = -14.4823°) or (b = 17.689814°,/ = 
-2.039549°), obtained from the HST archive (HST 2002). We 
detected 3518 stars in both colors in the HST images correspond- 
ing to the EROS monitored field; we systematically tried to as- 
sociate these stars with an EROS object within 1 arcsec. All of 
the 869 EROS -objects properly identified as stars in the field 
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Fig. 1. The Galactic plane fields (Galactic coordinates) monitored by EROS superimposed on the image of the Milky- way. The locations of our 
fields towards the spiral arms, as well as our Galactic bulge fields (not discussed in this paper) are shown. The large blue dot towards y Set indicates 
the position of the HST field used to estimate our star detection efficiency (see text). North is up, East is left. 




Fig. 3. (a) The Reros composite image (used to detect the cataloged 
stars) and (b) the U6FQ1104B-HST image of the same sub-field to- 
wards gs201. 



were associated with one or more HST star, allowing a study of 
the blending and of the detection efficiency as a function of the 
magnitude. We found that 56% of the EROS objects are blends 
with more than one HST-star within 1 arcsec distance. In this 
case the brightest HST-star accounts for an average of 72% of 
the Beros flux of the EROS object. These numbers vary with 
the EROS object magnitude as follows: 



Beros 


15-17 


17-19 


19-21 


fraction of blended EROS-objects 


70% 


59% 


54% 


contribution of main HST star 


88% 


11% 


68% 



The comparison of these numbers with the ones found from 
a similar study of a LMC dense field (Tisserand et al 2007) in- 
dicates that EROS cataloged objects towards gs201 are on av- 
erage less blended than the objects found towards dense re- 
gions of LMC. Such blending could affect the microlensing op- 
tical depth determination, as discussed in detail for SMC fields 
in (Afonso et al 2003), and the distribution of the lensing time 
scale tE (Rahvar 2004; Bennett 2005). The average densities of 
the EROS catalogs are similar towards the SMC and the spiral 
arms; therefore one should expect diflTerences in blending only if 
there is a difference between the spatial repartitions of the 2 stel- 



lar populations. Since our studies of HST images have shown 
that the spiral arm stars are less blended than the LMC stars, 
and considering the similarity between the SMC and the LMC 
populations, we conclude that blending should have less impact 
towards the spiral arms than towards the SMC. Therefore, to be 
conservative, we will use the estimates of (Afonso et al. 2003) 
as upper limits on the optical depth systematic uncertainties in 
Section 7.2. 

From the HST-EROS star association, we have extracted our 
detection efficiencies as a function of the Beros stellar magni- 
tudes (see figure 4). As F814W and F606W HST-WFPC2 filters 
are respectively very close to our Reros and Beros bands, we 
could directly measure detection efficiencies for HST objects. 
We found that every HST star that is detected in the EROS im- 
ages (i.e. that is located within 1 arcsec of an EROS object) in 
the Beros band is automatically detected in the Reros band (the 
reverse is false). This is due to the different limit magnitudes 
of the Beros and Reros templates. Therefore the efficiency to 
detect a HST star in EROS is the probability for that star to be 
found in the Beros band. The color-magnitude diagrams of Fig. 
5 show that the diagonal delimitations of the populations in the 
bottom right sector follow a Beros = constant line, thus con- 
firming that the detection threshold is set by Beros • We estimate 
the efficiency within the active region of the CCD-array, corre- 
sponding to the effective field of 0.71 deg^ for the full mosaic. 
We provide in Fig. 4 the probability for a HST star to be the main 
contributor of an EROS object. A star can also have a minor con- 
tribution to the flux of an EROS object, as a result of blending; 
we show also the probability for HST stars to contribute to an 
EROS object (even if not as the main contributor). 

4.2. The color-magnitude diagram 

Figure 5 gives the color-magnitude diagrams nerosiL V - /) of 
our catalogs^. The global pattern of these diagrams follows the 
expected magnitude versus color lines resulting from the light 
absorption of a distance-distributed stellar population. Two par- 
allel features are visible, with very different densities. 

^ 2D-tables of these diagrams can be found on the Web-site: 
http ://user s . lal . in2p3 . fr/moniez/ 
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Fig. 5. Color-magnitude diagrams UerosiL V - I) of our catalogs towards the 4 monitored directions. The grey scale gives the number density of 
stars per square degree, per unit of magnitude and per unit of color index. 



We were able to qualitatively reproduce these features with 
a (simple) simulated catalog (Fig. 21). The color-magnitude dia- 
gram of this synthesized catalog shows two parallel features due 
to the main sequence and the red giant clump, that are similar 
to the ones observed in the data. Without spectroscopic data or 
a more detailed simulation, it is not possible to go further than 
this qualitative comparison for the interpretation of the observed 
color-magnitude diagrams. 



4.3. Photometric precision 

To complete the description of our observations. Fig. 6 gives the 
average point-to-point photometric dispersion along the light- 
curves as a function of the magnitude /c. 



5. The search for lensed stars 

Our microlensing event detection scheme is the same as the one 
described in papers I and II. In the following, we will outline 
the few specificities that arise because of analysis improvements, 
specific seasonal conditions or particular problems, and because 
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HST 



i2 10 - 



EROS 



'EROS 



21 20 19 18 17 16 15 




21 20 19 18 17 16 15 



Fig. 4. - Top panel: The EROS (thick line) and HST (thin line) Beros = 
Ic + 0.6(^7 - Ic) magnitude distributions of the identified objects in a 
sub-field of gs201. Objects brighter than Beros = 16 are all identified in 
both images, but their magnitudes are systematically overestimated by 
our photometry in the HST image, explaining the apparent deficiency 
of bright HST objects. 

- Lower panel: The thin line shows the probability for an HST star to 
contribute to an EROS object, i.e. to be closer than 1 arcsec from such 
an object versus Beros - Ic + 0.6(^7 - Ic). 

The thick line gives the probability for an HST star to be the main con- 
tributor to the flux of an EROS object found within 1 arcsec. 



of the fact that the time baseline is twice to three times longer 
than in our previous publications. 




Fig. 6. Average photometric point-to-point precision along the light- 
curves versus Ic. The vertical bars show the dispersion of this precision 
in our source sample. The histogram shows the magnitude distribution 
of the full catalog (average over 4 directions). 



5.1. Prefiltering 

We used the same non specific prefiltering described in paper II, 
and preselected the most variable light curves satisfying at least 
one of the following criteria: 

- The strongest fluctuation along the light curve (a series of 
consecutive flux measurements that lie below or above the 
"base flux", i.e. the average flux calculated in time regions 
devoid of significant fluctuations) has a small probability 
(typically smaller than 10"^^) to happen for a stable star, as- 
suming Gaussian errors; 

- The dispersion of the flux measurements is significantly 
larger than expected from the photometric precision; 

- The distribution of the deviations with respect to the base 
flux is incompatible with the distribution expected from the 
measurements of a stable source with Gaussian errors (using 
the Kolmogorov-Smirnov test). 

The thresholds of these three criteria have been tuned to se- 
lect a total of ~ 20% of the light curves. After this prefiltering, 
2446843 light curves are entering the more discriminating anal- 
ysis described below. We also included a randomly selected set 
of light curves (~ 2%) to produce unbiased color-magnitude dia- 
grams and for our efiftciency calculation (Sect. 7.1). Furthermore, 
we have corrected the photometric measurements presenting a 
significant correlation between the flux and the seeing in a way 
that is described in (Tisserand 2004). 

5.2. Filtering 

- As in paper II, we first searched for bumps in each light 
curve. A bump is defined as a series of consecutive flux mea- 
surements that starts with a positive fluctuation of more than 
one standard deviation (+lcr) from the base flux, ends when 
3 consecutive measurements lie below Icr from the base flux 
and contains at least four measurements deviating by more 
than -\-lcr. We characterize such a bump by the parameter 
Q = -logio(P) where P is the probability that the bump be 
due to an accidental occurrence in a stable star light curve, 
assuming Gaussian errors. We select the light curves whose 
most significant fluctuation (bump 1) is positive in both col- 
ors. 

- Then we require the time overlap between the main bumps in 
each color to be at least 10% of the combined time intervals 
of the two bumps. 

- To reject most of the periodic or irregular variable stars, we 
remove those light curves that have a second bump (bump 2, 
positive or negative) with Q2 > Qi/2 in one color. 

After this filtering, the 1097 remaining light curves can 
be fitted assuming the simplest microlensing hypothesis, i.e. sl 
point-like source and a point-like deflector with a constant speed. 

5.3. Candidate selection 

The observed flux versus time data ^obsiU) is fitted with the ex- 
pression 0(0 = <^base X ^(0. where ^tase is the unmagnified 
flux and A(t) is given by expression (3). The candidate selection 
is based on the fit quality (x^) and on variables obtained from 
the ^base, to, tE and wq fitted parameters. We apply the follow- 
ing criteria, tuned to select not only the "simple" microlensing 
events, but also events that are aflTected by small deviations due 
to parallax, source extension, binary lens eflTects... mentioned in 
Sect. 2. The efiftciency to detect caustics should be very limited 
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with this set of cuts, but none was found from a systematic visual 
inspection of the 1097 light curves. 

- CI. Minimum observation of the unmagnified epoch : 

We first reduce the background due to instrumental eff'ects 
and to field crowding problems by selecting light curves that 
are sufficiently sampled both during the unmagnified and the 
magnified stages. For this purpose we define the "high" mag- 
nification epoch (called peak, labeled "t/ < 2") as the period 
of time during which the fitted magnification A is above 1 .06, 
associated to an impact parameter u < 2. The complemen- 
tary "low" magnification epochs, during which A < 1.06, are 
labeled ''base'' . We require that 



Mobs - ^Tu<2 > 600 days. 



(10) 



where Mobs = 2325. days is the observation duration, and 
ATu<2 is the duration of the "high" magnification epoch. 

C2. Sampling during the magnified epoch : We also re- 
quire that the interval between the peak magnification time 
^0 and the nearest measurement is smaller than 0.4 x ATu<2- 
C3. Goodness of a simple microlensing fit : To ensure the 
fit quality, we require x mi/ ^dof < 1-8 separately for both 
colors, where and the number of degrees of freedom Ndof 
are obtained from the full light-curve. 
C4. Impact parameter : We also require that the fitted im- 
pact parameter uq be less than 1 for both colors. 
C5. Stability of the unmagnified object : One important 
feature of a microlensing light curve is its stability during the 
low magnification epochs, except for the rare configurations 
of microlensed variable stars. We reject light curves with 



4ase(^)^xLe(B) 
Ndof{R)+Ndof{B) 



(11) 



where the xlase ^dof values correspond to the measure- 
ments obtained during the low magnification epochs. 
C6. Improvement brought by the microlensing fit com- 
pared to a constant fit : We use the same l^x^ variables 
as in paper II to select light curves for which a simple mi- 
crolensing fit is significantly better than a constant value fit: 



A 2 '^cst 



1 



(12) 



We select light curves with A;^| + A;^^ > 60. 
C7. Overlap in the two colors : Defining /STu<\ as the time 
interval during which the fitted magnification is larger than 
1.34 (t/ < 1), we require a minimum overlap between the 
time intervals found in the two colors: 



Mu<x{R)C\Mu<x{B) 
ATu<muATu<i(B) 



>0.4. 



(13) 



This loose requirement on the simultaneity of the magnifica- 
tions in the two colors allows one to keep a good sensitivity 
to "complex" microlensing events; for example, this cut tol- 
erates some diff'erence between the fitted impact parameters 
obtained in the two colors (which may occur in the case of 
strong blending). 

The number of microlensing candidates so far is 27 including an 
uncertain one, labeled GSA-ul (see below). The Ic magnitudes 
and (Vj - Ic) colors versus uq of these candidates are shown in 



Fig. 7 together with a sample of points representing the popula- 
tion obtained after selection of simulated events as explained in 
Sect. 7.1. One clearly sees how the maximal source magnitude 
required for detection decreases when the impact parameter in- 
creases. Annex A shows the light-curves and the finding charts 
of the 27 candidates, and table 2 gives their characteristics. The 
finding charts are obtained from the reference images used for 
the production of the catalogs. 

5.4. Non standard microlensing events 

Some of our candidates are significantly better fitted with mi- 
crolensing curves resulting from complex configurations than 
with the basic point-like source, point-like deflector with a 
constant- speed microlensing curve. The refinements that have 
been introduced in these cases are: 

- The blending of the lensed source with a nearby, unresolved 
object. In that case, the fight-curve (i>obs(ti) has to be fitted by 
the following expression 



0(0 = CX ^base X A(t) (1 - C) X ^ba 



(14) 



where C depends on the color. In supplement to the stan- 
dard fit, this fit provides the Cr and Cb parameters, where 
C=(base flux of magnified component)/(total base flux). In 
the notes of table 2, we give the magnitudes and colors of 
the microlensed components that take into account the color 
equations (9). 

- Parallax. Due to the rotation of the Earth around the Sun, 
the apparent trajectory of the deflector with respect to the 
line of sight is a cycloid instead of a straight line. For 
some configurations (a nearby deflector and an event that 
lasts a few months), the resulting magnification versus time 
curve may be affected by this parallax eff'ect (Gould 1992, 
Hardy & Walker 1995). The specific parameters that can be 
fitted in this case are the Einstein radius and an orienta- 
tion angle, both projected on the observer's plane which is 
orthogonal to the line of sight. 

- "Xallarap". This eff'ect is due to the rotation of the source 
around the center-of-mass of a multiple system. In this 
case, the light-curve exhibits modulations with a charac- 
teristic time given by the period of the source rotation 
(DeruQetaL 1999, MoUerach & Roulet 2002). Assuming a 
circular orbit, the extra-parameters to be fitted or estimated 
are the orbital period P, the luminosity ratio of the lensed ob- 
ject to the multiple system, and the projected orbit radius in 
the deflector's plane p = ax/ Re, where a is the orbit radius 
and X = Diens/D 

source • 

In some cases, the values obtained for tE with the basic fit and the 
refined one may differ considerably. For time duration studies 
and for optical depth calculations, we use the tE values given 
by the best fit. But as far as efficiency values are concerned, we 
must use those obtained with the standard fit since they are the 
ones that enter the selection procedure. 



6. The microlensing candidates 

6.1. General features 

In order to quantify the relevance of the interpretation of the 27 
selected objects as microlensing events, we define two variables 
as follows: 
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Table 2. Characteristics of the 27 microlensing candidates. For those events that have a better fit than the point-like point-source constant speed 
microlensing fit (the so-called standard fit), we also provide the standard fit parameters. 

- Names in bold type correspond to events selected for the optical depth and duration analysis (with uq < 0.7). 

- Ic (Vj - Ic) are the fitted unmagnified magnitudes of the lensed object (including the contribution of a possible blend). 

- ^0 is the time of maximum magnification, given in HJD-2,450,000. 

- tE is the Einstein disk crossing time, in days. 

- uq is the dimensionless impact parameter. 
-X^/dof corresponds to the best microlensing fit. 

- T is the individual contribution of each event to the optical depth towards the corresponding target. In the case of "non standard" events, we use 
the tE value obtained from the best (non standard) fit and the efficiency evaluated at tE of the standard fit (see text) for the calculation of r. 



candidate 


field 




6° (J2000) 


Ic(Vj-Ic) 


fo(days) 


^^(days) 


Uo 


x'/d.o.f 


t(10-^) 


note 


7 Set 


GSAl 


200 


277.2888 


-14.2528 


18.3 (3.1) 


301.2+ 0.1 


64.0+ 1.2 


.043+.0010 


299.7/435 


0.146 




GSA8 


200 


276.8042 


-15.0311 


16.6 (3.7) 


996.9+ 0.1 


40.6+ 1.0 


.145+.003 


981./533 


0.114 


(1) 








Standard f 


it parameters: 


993.0+ 0.1 


35.2+ 0.7 


.155+.002 


1400./535 






GSA9 


200 


277.1750 


-15.1644 


18.8 (2.8) 


1760.3+ 1.7 


57.9+ 3.6 


.482+.0158 


262.9/565 


0.142 




GSAIO 


200 


277.2813 


-14.8931 


18.1 (2.3) 


1806.1+ 0.8 


24.6+ 1.3 


.574+.0179 


237.9/596 


0.081 




GSAll 


201 


278.1650 


-14.1094 


18.8 (2.5) 


1725.3+ 0.4 


44.3+ 1.4 


.187+.0049 


367.2/558 


0.116 




GSA12 


203 


278.5875 


-13.9794 


17.1 (1.8) 


1378.6+ 0.2 


50.1+0.7 


.225+.0030 


89.3/359 


0.123 




GSA13 


203 


278.9404 


-14.5803 


17.2(1.9) 


313.9+ 1.2 


37.2+ 2.1 


.898+.0153 


244.9/617 






GSA14 


204 


278.4388 


-12.8678 


16.7 (2.2) 


1637.7+ 3.4 


68.4+ 3.7 


.785+.0097 


421.2/392 




(2) 




GSA15 


301 


281.0654 


-6.0339 


18.3 (3.5) 


1399.8+ 1.4 


72.2+ 2.8 


.337+.0126 


212.2/411 


0.110 


edge 


GSA16 


301 


280.7646 


-6.7583 


16.4 (3.1) 


1997.0+ 3.2 


60.6+ 4.0 


.796+.0141 


341.8/361 






GSA17 


302 


281.3950 


-7.8867 


16.4(1.7) 


1947.2+ 3.8 


50.0+ 2.4 


.532+.1079 


156.9/400 


0.096 




GSA18 


304 


282.2879 


-7.2500 


16.1 (2.1) 


1718.7+ 0.1 


55.0+ 2.0 


.137+.0009 


133./514 


0.098 


(3) 








Standard f 


it parameters: 


1718.4+ 0.1 


58.0+ 0.3 


.137+.0009 


155.6/516 






7 Nor 


GSA2 


400 


242.9592 


-52.9464 


18.6 (2.3) 


534.4+ 0.2 


98.3+ 0.9 


.342+.002 


973.4/934 


0.059 


(4) 








Standard f 


it parameters: 


533.6+ 0.5 


137.8+ 2.6 


.233+.0029 


1196.5/937 






GSA19 


401 


244.1379 


-52.0272 


15.5 (5.1) 


2367.7+ 1.3 


90.4+ 3.0 


.043+.025 


893./880 


0.063 


(5) 








Standard f 


it parameters: 


2373.5+ 0.1 


93.1+0.7 


.022+.007 


1388.5/888 






GSA20 


402 


243.7758 


-52.9700 


16.4(1.3) 


2465.5+ 1.0 


40.+ 5.0 


.72+.02 


414./696 


0.039 


(6) 








Standard f 


it parameters: 


2487.1+ 0.3 


46.3+ 0.7 


.565+.0050 


712.4/698 






GSA21 


404 


244.3063 


-53.1100 


16.8 (2.7) 


1587.3+ .03 


74.+ 3.0 


.0142+.0008 


259./565 


0.077 


(7) 








Standard f 


it parameters: 


1587.2+ .03 


39.1+0.2 


.037+.0007 


1884.4/567 






GSA22 


404 


244.4263 


-54.0508 


18.4 (2.0) 


2182.4+ 0.2 


26.6+ 1.1 


.048+.0180 


429.5/742 


0.031 




GSA23 


404 


245.1208 


-53.9825 


18.3 (1.7) 


1573.8+ 3.8 


78.5+ 5.7 


.542+.0152 


522.3/784 


0.062 


corner 


GSA24 


406 


246.5442 


-54.0394 


18.0(1.9) 


2002.5+ 1.4 


55.5+ 2.6 


.720+.0158 


579.4/786 






GSA25 


408 


247.6917 


-53.9281 


21.1 (1.4) 


850.9+ .03 


67.6+ 2.9 


.003+.0001 


876.2/771 


0.057 


(8) 


GSA3 


409 


244.1129 


-54.6303 


17.7 (1.4) 


696.0+ 2.0 


60.4+ 3.0 


.615+.0102 


606.7/1090 


0.051 




GSA26 


411 


241.8729 


-55.3814 


17.8 (2.2) 


1642.1+ 0.3 


23.2+ 0.8 


.504+.0138 


441.5/759 


0.030 




GSA27 


411 


242.4846 


-55.2292 


18.3 (1.7) 


2193.8+ 0.1 


6.8+ 0.4 


.210+.0068 


433.6/831 


0.022 




OWlus 


GSA28 


501 


202.2838 


-64.2750 


19.3 (3.7) 


1992.2+ 0.4 


205.+ 20.0 


.029+.004 


717/499 


0.431 


(9) 








Standard fit parameters: 


1992.0+ 0.4 


87.3+ 3.0 


.094+.0046 


868.5/500 






GSA29 


502 


204.0683 


-63.7117 


19.3 (2.5) 


1229.7+ 0.3 


74.2+ 2.7 


.082+.0042 


161.7/354 


0.166 




GSA30 


505 


199.2942 


-64.2592 


16.5 (2.7) 


2396.9+ 0.1 


12.4+ 0.2 


.062+.0023 


792.0/856 


0.073 




Uncertain eandidate 


GSAul 


202 


278.0371 


-13.2851 


18.1 (2.9) 


1695.8+ 6.9 


409.3+20.9 


.708+.0155 


426.9/613 







Notes: (1) GSA8: Blended; Cr = 1.00 + 0.03, Cb = 0.68 + 0.02; lensed star has /* (V; - /*) = 16.6 (4.4). (2) GSA14: Light-curve exhibits 
typical features of a binary lens system. Given the small number of measurements with significant magnification, no reliable analysis of the shape 
can be performed. (3) GSA18: Parallax; projected Einstein radius in the solar plane - 12.5 + 7.0AU. (4) GSA2: Described in paper I. Found at 
that time as the first candidate for a binary lensed source (Xallarap). (5) GSA19: Xallarap and blend; the best fit is performed ignoring the 3 most 
magnified measurements, that are afi'ected by the non-linearity of the CCD. Cr = I., Cb = 0.160 + 0.013; lensed star has /* (V* - /*) = 15.5 (8.4). 
The light-curve distortion could be due to the face-on circular orbiting of the source around the center-of-mass of a system including a non 
luminous object, with period Pq = 294. + 47. days, and with a projected orbit radius of p = ax/ Re = 0.081 + 0.023, where a is the orbit radius 
and X = Diens ID source- Scc also text. (6) GSA20: Parallax; ve = 0.94 + 0.07At/. (7) GSA21: Blended; Cr = 0.53 + 0.02, Cb = 0.34 + 0.02; lensed 
star has /J (V} - /J) = 17.5 (3.5). (8) GSA25: An improbable configuration, but a genuine one (very small uq on a very faint star). (9) GSA28: 
Blended; Cr = 0.30 + 0.03, Cb = 1; lensed star has /* (V} - F^) = 20.6 (1.5). Thox^/dof of the fit is aff'ected by an underestimate of the errors 
due to bright neighboring stars. 



- Ideally, the goodness of the microlensing fit should be uni- 
form throughout the observation duration. Here we use fits 
made separately in the two colors. Let ^^^2 ^u<2 be the 
microlensing fit and the number of degrees of freedom, 



restricted to the high magnification epoch (u < 2, A > 1.06, 
see Sect. 5.3). ^^^xlase ^^ase be the complementary vari- 
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Fig. 7. Top panel: Ic versus fitted uq for the 27 microlensing candidates. 
Bottom panel: Vj - Ic versus fitted uq. uq is the fitted value assuming a 
point-like source and a point-like deflector with a constant speed. The 
red dots (respectively black dots) correspond to events with uq < 0.7 
used for optical depth studies (resp. 0.7 < uq < 1.). The small dots are 
the simulated events that satisfy the microlensing selection criteria. 



ables outside the peak period. The variable 



nu<2{R) + nu<2(B) maseiR) + nbase(B) 



Fig. 8. 6 fit versus (A;^^ + Ax^)/tE for the microlensing candidates. Red 
dots correspond to events with fitted uq < 0.7. The small dots represent 
the simulated events that satisfy the microlensing selection. 



1 



1 



nuKliR) + nu<2(B) flbaseiR) + nbase(B) 



(15) 



quantifies the diflference of the standard fit quality during and 
outside the microlensing peak, expressed in standard devia- 
tions (thanks to the second factor). A negative value of Sftt 
(< -5) is an indication of a non constant base, and points 
to a variable star instead of a microlensing event. For non- 
standard microlensing (parallax, blending...) S fit will be pos- 
itive and may be large (> 10), because the fit is expected to 
be less good in the peak than in the base. 
Many of the EROS instrumental defects — such as bad pix- 
els or diflTraction features — have a long lifetime, and last for 
entire observing seasons. This produces long time scale false 
candidates. A signal to noise indicator is provided by the ra- 
tio (A;^| -h Axj^)/(tE/^ day) where Ax^ is defined above (cri- 
terion C6) and where tE characterizes the event time scale. 



Figure 8 shows the distribution of Sfit versus (Ax^ + 
Axj^)/(tE/l day) for the data satisfying the filtering conditions, 
for the final candidates and for the simulated sample (see Sect. 
7.1). Within our final sample, three subgroups are apparent: 

- a few (6) selected events have both a large positive S fit and 
(A;^^ -h Ax^)/tE. These are events for which a non-standard 
microlensing fit provides a better interpretation. Each of 
them is discussed in the remarks of table 2. 

- the bulk of our final sample (20) are events with a large 
(Ax^Q + Axj^)/tE and Sfit compatible with 0, as expected for 
standard microlensing events (and as is the case for our sim- 
ulated sample). 

- Event GSA-ul has a small (A;^^ -h Axj^)/tE. After visual in- 
spection (see Annex A, last event), we cannot exclude a mi- 
crolensing interpretation, but the long duration and the lack 
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of a reliable base make it very uncertain, considering the rel- 
atively low value of A;^^ + A;^^ (only 73.). The status of this 
candidate remains pending until further observations over a 
longer time range can be made. One should keep in mind 
that confirmed events of this type would give a major con- 
tribution to the optical depth (GSAul would contribute for 
~ 0.5 X 10"^ towards 6 Mus). 

6.2. Comparison with the EROS 3 year analysis (paper II) 

We first checked the coherence between the present results and 
those of paper II. Three additional candidates (GSA8, 13 and 25) 
with a maximum occurring during the first three years have been 
found. GSA13 and GSA25 belong to subfields that were not ana- 
lyzed in paper II. GSA8 is located at the border of two subfields, 
and was missed by our previous analysis that did not systemati- 
cally explore the overlapping regions between subfields. 

Four of the 7 candidates, all towards y Set, found in paper II 
are now rejected for the following reasons: 

- GSA4 and GSA7 both showed a second fluctuation after the 
first three years. 

- For GSA5, the;^^ improvement when replacing a constant fit 
by a microlensing fit is no longer significant enough, due to 
the low signal to noise ratio that prevails for its light curve 
during the 7 years of data taking. 

- GSA6 was found to have an impact parameter of 0.98 ± 0.04 
in paper II. Taking into account the full light curve, the new 
fitted value is uq = 1.03 ± 0.07, now just above our thresh- 
old. Incidentally, Ax^ is also much smaller than our threshold 
(60), indicating that the previous selection of this event could 
have been due to a fluctuation. 

One notices that these rejected candidates were the low sig- 
nal/noise ones towards y Set. Clearly, 7 years of observations 
allow a much better noise reduction than 3 years. 

6.3. Overlap with other published surveys 

A very small region of y Nor overlaps with the OGLE II mi- 
crolensing survey (OGLE webpage). No event from this region 
was reported in the latter survey (Udalski et ah 2000a). 

A small region of our survey overlaps the MACHO fields 
(Thomas et al 2005). Amongst the 9 MACHO candidates or 
alerts found around y Set, 3 are located within one of our mon- 
itored fields, but have not been selected in our analysis for the 
following reasons: 

- MACHO alert number 302.44928.3523 is too faint to be 
measured in Beros and no measurement was made in Reros 
within 40 days of the magnification maximum. Nevertheless, 
an object clearly appears in the Beros images around the 
maximum magnification date. 

- MACHO alert number 301.45445.840 is too faint to be in 
the EROS catalog. Furthermore, EROS missed the event as 
its time of maximum magnification was 106 days before the 
first EROS observation of the corresponding field. 

- MACHO alert number 302.45258.1038 was very close to 
one of the gaps located between CCDs. Thus many measure- 
ments are missing. Our standard procedure does not try to 
recover complete light-curves in such a case, and the stan- 
dard light-curve failed our selection process. Nevertheless, 
we confirm the presence of the bump at the right time, with 
the maximum magnification occurring during the very first 
days of the EROS data taking. 
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Fig. 9. to versus fitted uq for the simulated events satisfying the analysis 
criteria (small dots) and for the detected candidates (big dots). Red dots 
correspond to events with fitted uq < 0.7 ("standard" fit). In the case of 
complex events, the best fit uq value is plotted. 



6.4. Statistical properties of the candidate parameters 

6.4.1. The lens configurations 

Microlensing events occur with a flat-distributed impact param- 
eter and minimum approach time. The sample of observed mi- 
crolensing event (^o, ^^o) configurations should be statistically 
representative of such a distribution after taking into account our 
detection eflftciencies. This is illustrated in Fig. 9, where simu- 
lated events are generated as described in Sect. 7.1. 

6.4.2. The microlensed star population 

The microlensed star population should also be representative 
of the monitored population weighted by the microlensing de- 
tection efficiencies and by the optical depth that may vary from 
source to source. As the sources are likely to be distributed along 
the line of sight, a possible variation of the optical depth with 
distance must be considered in the data analysis. As the light of 
a remote source is expected to be more reddened than the light 
of a close one, the optical depth r should increase on average 
with the color index. Figure 10 shows our color-magnitude di- 
agram, weighted by the microlensing efficiencies and assuming 
the same optical depth for all stars. It is directly obtained from 
the simulated events that satisfy the analysis requirements. The 
distribution of the observed candidates is less peaked than the 
simulated one in the low color index region because the most 
reddened stars are more likely to be lensed. We were able to 
qualitatively confirm this color bias through the catalog pro- 
duced with a simple simulation towards y Set described in Sect. 
9.2, that takes into account the source distance distribution (Fig. 
11 left). Fig. 11 (right) shows that the color distribution of the 
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Fig. 11. Left: Optical depth as a function of the distance (thick line), 
source distance distribution of a simulated catalog (thick histogram), 
and distance distribution weighted by the optical depth (thin histogram). 
Right: color distribution of the stars of the simulated catalog (thin line) 
and expected color distribution of the lensed stars (thick line). 



tation will be provided in a forthcoming publication accordingly 
to the guidelines given in Sect. 9. 

Two outliers need a specific comment. On closer inspection, 
it appears that GSA25 is a genuine microlensing candidate of a 
very faint star. It was detected because of the very strong mag- 
nification. This is a rare case, but there is no reason to discard it 
from our list. GSA19 is a very bright and very red object. It could 
be a strongly absorbed nearby star lensed by a closer object. 

The spatial distribution of the candidates shown in Fig. 12 
does not indicate any remarkable concentration. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 
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Fig. 10. Color-magnitude diagram and projections of the simulated 
events satisfying the analysis criteria (small dots) and the detected can- 
didates (big dots). The arrows show {Vj - Ic, Ic) of the magnified com- 
ponent in the case of blending (see notes of table 2). The red dots rep- 
resent those events that are used for the optical depth estimates. The 
histograms of these events are superimposed on the projections (not 
normalized). 



lensed sources (obtained by weighting with the optical depth) 
is significantly biased towards the red color with respect to the 
simulated distribution of detected sources. We conclude that the 
distance scattering of the sources can explain the observed bias 
of the lensed stars towards red color. A more complete interpre- 



6.5. Domain of sensitivity of the analysis 

The CI and C2 cuts that use the ATu<2 duration of the stage with 
magnification A > 1.06 mainly aff'ect the light-curves which 
show long bumps. Removing these two cuts adds 6 candidates, 
all of very long duration, that have a poor signal to noise ratio. 
As is the case for the GSA-ul event, only a very long monitor- 
ing could change the status of such candidates and improve or 
degrade their signal to noise ratio. Therefore one should keep 
in mind that the optical depths we publish in this paper are al- 
most insensitive to events with tE > 700 days (cf. the detection 
efficiency versus tE curve in Fig. 13). 

7. Optical depth 

To obtain reliable optical depth values, we use a sub- sample 
of good quality candidates which should be almost free 
of the microlensing like variable objects that have been 
identified towards other EROS targets (Tisserand et al 2007; 
Hamadache et al 2006) and in Sect. 6.2. For this purpose, we 
will only keep those candidates that have uq < 0.7 in the fit 
that assumes a point-like source and a point-like deflector with 
a constant speed. This is approximately equivalent to requesting 
A^ax > 1.68. Events GSA13, 14, 16, 24 and GSAul are then 
discarded for the optical depth analysis. 

7.1. Microlensing detection efficiency 

We present here the efficiency calculation for the detection of 
events with uq < 0.7. As for our previous papers, we calculate 
our detection efficiency by superimposing simulated events on 
measured light curves from an unbiased sub- sample of our cat- 
alog. Events are simulated as point-source, point-lens constant 
velocity microlensing events, with parameters uniformly span- 
ning a domain largely exceeding the domain of EROS sensitivity 
(uq up to 2, 1 day < tE < 900 days, to generated from 150 days 
before the first observation to 150 days after the last). Efficiency 
is defined as the ratio of events satisfying the selection cuts to the 
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Fig. 12. From left to right: spatial distribution (a, 6) of the 27 candidates in the monitored fields towards Mus, y Nor, y Set and j3 Set. East is 
right, North is up. The area of a marker is proportional to the contribution of the candidate to the optical depth towards the corresponding target. 
Events with uq > 0.7 are marked with a star. 



number of events generated up to uq = 1. Figure 13 (upper left 
panel) shows the EROS efficiency as a function of the source 
position in the color-magnitude diagram averaged over all the 
other parameters and over all directions. The other frames of 
Fig. 13 show the efficiency as a 2D-f unction of Ic of the lensed 
star and tE, and as a ID-function of tE,uo, Ic and Vj - Ic, aver- 
aged over all the other parameters, for each monitored direction. 
The efficiency is significantly better towards y Nor, because of 
the higher sampling of the light curves in this direction. 

7.2. Optical depth determination 

The optical depth values obtained from the 22 events that satisfy 
Uq < 0.7 are given in the first part of table 3. The average over 
all directions is defined as the proportion of stars covered by an 
Einstein disk. It is given by 



T fields = - With n = - 



1 



tE 



2N[M\ e\tE) 

* obs events 



where A^^, AP^^, e\tE) are respectively the number of stars 
monitored, the observation duration and the microlensing selec- 
tion efficiency towards direction /. As usual, selection efficiency 
is relative to events with U{) < I (even though the efficiency is 
almost zero for uo > 0.7). As explained at the end of Sect. 5.4, 
we use ^(j^lsumdard\it)) cxprcssion of T for non-standard 

microlensing events. The statistical uncertainties are estimated 
from the definition of the classical 68% confidence intervals 
(Feldman & Cousins 1998), multiplied by the factor 



^< tlle\tE) > 
< tEle{tE) > 



following Han & Gould 1995. According to the discussion of 
Sect. 4.1, we assume a 10% maximum systematic uncertainty 
due to blending eflTects, as in Afonso et al 2003, and we account 
for another 5% uncertainty due to the statistical limitations in 
the determination of the efficiencies. These errors are small com- 
pared to the smallest statistical uncertainty (~ 30%) estimated in 
this paper. 





6* Mus 


7 Nor 


7 Set 








-1.46 


-2.42 


-2.09 


-2.15 




1° 


306.56 


331.09 


18.51 


26.60 




Observations 


All 


TX 10^ 


•'^^-.52 


49+21 


79+.41 
• ^^-.28 


••^ -.20 


51+-^^ 

•^^-.13 


events 


3 


10 


6 


3 


22 


fE (days) 


97 + 47 


57 + 10 


47 + 6 


59 + 6 


60 + 9 




80 


29 


13 


10 


40 


6 corrected tE 


65 + 45 


43 + 10 


45 + 6 


58 + 6 


48 + 9 


€ corrected cr^^ 


75 


31 


13 


9 


38 


median tE 


74.2 


64. 


47.2 


55. 


56.5 



T (xlO^) from models 




A 


0.32 


0.48 


0.79 


0.60 


2.3 


+ spiral 


0.56 


0.69 


1.07 


0.83 


7.9 


B 


0.34 


0.51 


0.85 


0.64 


2.9 


+ spiral 


0.61 


0.72 


1.13 


0.90 


10.2 


C 


0.47 


0.78 


1.11 


0.95 


12.5 


+ spiral 


0.71 


1.18 


1.43 


1.23 


36.1 


D 


0.32 


0.56 


0.41 


0.38 


1.4 


model 1 


0.42 


0.52 


0.71 


0.57 


1.8 


model 2 


0.54 


0.68 


0.90 


0.74 


5.4 



Predictions from model 1 


^ events 
tE 

median tE 


2.8 9.9 7.1 6.3 
73.8 67.9 37.9 60.2 

63 54 36 48 
54.5 52.5 28.0 46.5 


tE from published models 


C (no effic.) 
C+ spiral 
D 


45 28 25 27 
38 38 44 40 
80.9 86.5 76.9 77.4 



Table 3. 

- Observed optical depth r, number of events N events, average tE, disper- 
(16) sion cr^g, efficiency corrected tE and cr^^ and median time scale for each 

monitored direction. Here, we consider only the events used for the op- 
tical depth estimates {i.e. with uq < 0.7). 

- Optical depth predictions from models A, B, C, including or not a 
spiral structure (see text), and from model D, and estimates up to 7 
kpc from the two simple Galactic models described in table 4. The;^^'s 
quantify the adequation with the observations (see Sect. 8). 

- Expected numbers of events and duration distribution parameters from 
model 1 and from models C and D. tE from model C do not take into 
account the detection efficiency. 



7.3. Robustness of the optical depth values 

We have studied the stability of the optical depth averaged over 
all fields by changing the selection cuts. Fig. 14 gives the vari- 
ation of f fields with the Ax^ threshold. Relaxing this threshold 



makes the optical depth increase rapidly, due to the inclusion 
of false events, as a logical result of the analysis optimization. 
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Fig. 13. Microlensing detection efficiency. 

Upper left: the efficiency in the (Vy - /c, /c) plane averaged over uo,to, Ie and over the 4 monitored fields. 
Upper right: the efficiency in the {Ie, Ic) plane averaged over wo, to, Vj - Ic and over the 4 monitored fields. 

Other panels: efficiency as a function of the lensed star magnitude Ic, as a function of the color index (Vj - Ic) and as a function of the event 
parameters tE and uq. Each curve gives the efficiency averaged over all the other parameters. 
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Fig. 14. Variation of the number of selected candidates (histogram, right 
scale) and of ffieUs (dots, left scale) with the Ax^ threshold. The red 
vertical line shows our cut. Only statistical errors are plotted. 



Fig. 16. Number of selected candidates (histogram, right scale) and 
"^fields (dots, left scale) for the sub-samples of stars brighter than the 
Ic threshold. Only statistical errors are plotted. 
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Fig. 15. Variation of the number of selected candidates (histogram, right 
scale) and of ffieUs (dots, left scale) with the uq threshold. The red verti- 
cal line corresponds to our selection. Only statistical errors are plotted. 



But using a stricter cut does not significantly change the optical 
depth values, as long as the statistics remains significant. 

Fig. 15 also gives f fields as a function of the uq threshold. 
Our result does not depend on this cut within statistical errors, 
showing that we probably found an optimum between the quality 
of the events and the number of those kept for our optical depth 
calculations. 



Fig. 16 and 17 show the variation of f fields with the maxi- 
mum magnitude 4 of the source population and with the min- 
imum color index Vj - Ic. There is no evidence for a vari- 
ation with Ic threshold. As discussed below in Sect. 9.3, this 
comes from the fact that the variation of the optical depth with 
the distance does not result in a variation with the magnitude 
as more distant identified sources do not appear fainter in our 
catalog on average. Interestingly one may use these figures to 
extract r for specific stellar populations, in particular the popu- 
lation of the brightest stars with Ic < 18.5 that are better identi- 
fied and measured, and that suff'er less blending (see Sect. 4.1). 
Moreover, our catalog is almost complete up to this magnitude, 
as our star detection efficiency is large (see Fig. 4). Using such a 
sub-catalog of bright stars should make the interpretation easier 
within a Galactic model framework as will be discussed here- 
after. 



8. Discussion: comparisons with simple models 

8.1. Optical depth 

We will consider here 4 published optical depth calculations and 
our own calculations based on simple Galactic models (model 
1 without a thick disk, model 2 with a thick disk) that we al- 
ready used for discussions in papers I and 11. The main revi- 
sion to these 2 models since our previous papers comes from 
the bulge inclination; we now take O = 45° instead of 15° 
(Hamadache et ah 2006 and Picaud & Robin 2004) as the angle 
of the outer bulge with respect to the line of sight towards the 
Galactic center. The only impact of this change is a little varia- 
tion of our optical depth value towards y Set. We also completely 
neglect any contribution from the halo to the optical depth, in the 
light of the latest EROS results towards the Magellanic Clouds 
(Tisserand et al. 2007). As in our previous papers, we performed 
simple optical depth calculations assuming all the sources to be 
at the same distance. More sophisticated modelling based on the 
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Fig. 17. Number of selected candidates (histogram, right scale) and 
fields (dots, left scale) for the sub-samples of stars redder than the color 
index (Vj - Ic) threshold. Only statistical errors are plotted. 



Table 4. Parameters of the Galactic models 1 and 2 used in this article. 





Parameter 


model 1 


model 2 




Ro (kpc) 


8.5 




E (Mopc-^) 


50 




H (kpc) 


0.325 




R (kpc) 


3.5 


Thin disk 




4.3 




(Tr {km s~^) 


34. 




0-0 (km s~^) 


28. 




(Tz (km s~^) 


20. 








35 




H (kpc) 




1.0 




R (kpc) 




3.5 


Thick disk 






3.1 




(Tr {km s~^) 




51. 




(Te {km s~^) 




38. 




(Tz {km s~^) 




35. 




a (kpc) 


1.49 




b (kpc) 


0.58 


Bulge 


c (kpc) 


0.40 




Inclination O 


45° 




Mb{xIO^^Mq) 


1.7 




o-buige {km ^"0 


110. 



Fig. 18. Expected optical depth up to 7 kpc at a Galactic latitude b = 
-2.5° for model 1 (thick curve) and model 2 (thin curve) which assumes 
an additional thick disk. The dotted lines show the excursion of the 
model 1 predictions when varying the latitude by +1°. The measured 
optical depths are given for our 4 targets. 



taken from Dwek et ah 1995 model G2, given in Cartesian coor- 
dinates: 



PB = 



G.Slnabc 



where Mb is the bulge mass, and a, b, c the length scale factors. 
Fig. 18 shows the measured optical depth as a function of the 
Galactic longitude, with the expectations from models 1 and 2 at 
7 kpc and at a Galactic latitude b = -2.5°. We also show in this 
figure the eff'ect of a Galactic latitude change by ±1°. 

The optical depth predictions of model 1 and 2 and of the 4 
following models are reported in table 3: 

- Model A, from Binney^^<2/. 1997, revised by 
Bissantz a/. 1997, has a cuspy and flat Galatic bar, 
inclinedby O ~ 20°. 

- Model B, described in Dwek et ah 1995, has a wider cuspy 
bar, inclined by O ~ 24°. 

- Model C, described in Freudenreich 1998, has a more ex- 
tended and difl'use bar, inclined by O ~ 14°. 



guidelines discussed in Section 9 will be considered in a forth- 
coming paper. We give in table 4 the list of the geometrical 
and kinematical parameters used in these models I and IL The 
disk densities are modeled by a double exponential expressed in 
cylindrical coordinates: 

m ^ ^ i-\z\ 
p.(7^,z)=-exp^^^jexp^— 

where X is the column density of the disk at the Sun position 
Rq, H the height scale and h the length scale of the disc. The 
density distribution for the bulge - a bar-like triaxial model - is 



The optical depths towards the directions monitored by EROS 
have been calculated for these 3 models and discussed by 
Evans & Belokurov 2002, who have also included the contribu- 
tion from the disk and considered separately the impact of a spi- 
ral structure. These 3 models have been normalized to have the 
same total mass (1.5 x lO^^M©) within 2.5 kpc. Therefore, they 
mainly difl'er by the shape details and the orientation of the bar. 

- Model D, described in Grenacher et al 1999, has a bar that 
is similar to the one of our models 1 or 2, but inclined by 
O ~ 20°, with a combination of a light thin disk plus a thick 
disk and a halo contribution. 
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The adequation of all these models with the observed optical 
depths T can be compared through the value of 



2 

model 



E 

targets i 



(Tiimodel) - Tiiobserved)) 



(17) 



where cr/ is the error interval of the t/ determination (as errors 
are asymmetrical, we consider the largest one for each measure- 
ment). The numbers reported in table 3 show that our model 1 
is clearly favored by the data, and also models A and D with- 
out spiral structure, which are the ones that predict the smallest 
optical depths. 

We cannot draw more conclusions about our model 1 , as we 
know that it is not a realistic description, since all targets are 
supposed to be at the same distance. One may also notice that 
the extrapolation of any bulge model to the relatively distant re- 
gion that we monitored is very uncertain. Nevertheless, it seems 
that "heavy" models trying to include a thick disk or any spiral 
structure are not favored. These results confirm the conclusions 
of Hamadache et al 2006, in particular for model C, that is also 
disfavored by the present data. 

8.2. Event duration distribution 

Fig. 19 gives the tE distributions for the observed events, com- 
pared to the expected distributions from model 1 . The procedure 
to build these distributions is the same as in paper II: 

- The mass function for the lenses is taken from 
Gould, Bahcall&Flynn 1997 for both the disk and the 
bulge. 

- The solar motion with respect to the disk is taken from 
(Delhaye 1965): 



VeR = -10.4, vo^ = 14.8, Vo, = 13 (km/s). 
- The global rotation of the disk is given by 



Vrotir) = Vrot,Q X 



1.00762 



0.0394 



-h 0.00712 



(18) 



(19) 



where V, 



rot,Q 



220 km/s (Brand & Blitz 1993). 



The peculiar velocity of disk stars is described by an 
anisotropic Gaussian distribution and a velocity dispersion 
given in table 4. 

The velocity distribution of the bulge stars is given by 



/r(vr) = 



1 



-Vt exp - 



bulge 

with (Thi^i^e ~ 110 km/s 



bulge ) 



(20) 



The Kolmogorov-Smirnov tests indicate that model 1 pro- 
duces Ie distributions that are fully compatible with the obser- 
vations, except for y Set where the approximation of a single 
distance for all the sources is the most questionable (due to the 
contribution of the sources belonging to the elongated bulge). 
Nevertheless, the supplement of short events towards this direc- 
tion that is visible in the Ie distribution corresponds to the expec- 
tations from bulge lenses. We found that our simulation results 
are not dramatically changed if we change the contribution of the 
lightest lens objects (around O.IM©) by 50%. This is mainly due 
to the low detection efficiency for the short events expected from 
low mass lenses (this loss of efficiency is specially noticiable for 
low mass bar lenses because their velocity is larger). 
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Fig. 19. measured Ie distributions compared to the one calculated with 
our model 1 for one day intervals, Pks gives the Kolmogorov-Smirnov 
test probability of compatibility of the shapes. Notice that the scales are 
different for Q Mus, The Ie values of the observed events are given by 
the vertical bars. 



Table 3 gives the observed averages, dispersions and medi- 
ans of Xe, with the predictions from model 1 and models C and D. 
Only mean duration predictions have been estimated for model 
C (with no efficiency corrections) and D (taking into account 
EROS microlensing detection efficiency). The general tendency 
that emerges from the observations and that is confirmed by the 
predictions is that increases with the Galactic longitude. This 
is due to the fact that the direction of global motion of the lenses 
tends to align with the line of sight when the longitude increases, 
and consequently the transverse speed becomes smaller. 

As the variation of the microlensing detection efficiency with 
Ie is not taken into account in the predictions of model C, we 
should compare here the predicted values with the efficiency 
corrected means given by 



Mle(tE) 



VeitE) 



(21) 



The observed corrected means are significantly larger than the 
predictions of model C with no spiral structure. This discrep- 
ancy could be explained (at least partially) by the absence of 
very short time scale events that have a negligible detection effi- 
ciency and can be totally missed in our statistically limited sam- 
ple; a more precise comparison could be done between the ob- 
served Xe distribution and a modified predicted distribution tak- 
ing into account the microlensing detection efficiency as a func- 
tion of Xe (Fig. 13). Keeping in mind the possibility of unde- 
tectable events, it seems that model C without spiral structure is 
disfavoured. Evans & Belokurov 2002 have emphasized the im- 
pact of streaming in the bar on the Xe distribution around the 
Galactic center and the impact of the spiral structure on X^ in ev- 
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Galactic longitude 1 (degree) 



Fig. 20. Model 1 predictions for the optical depth at 6, 7, 8, 9, 10, 1 1 and 
12 kpc (from lowest to highest curve) at a Galactic latitude b = -2.5° 
as a function of the Galactic longitude. The measured optical depths are 
given for our 4 targets. 



ery direction. Such a streaming and the spiral structure have a 
serious impact on the 4 values, and could reconcile the model 
expectations with the tE data, but lighter disk and spiral structure 
should be considered to fit the optical depth data as well. There 
is probably some margin of freedom for doing this, as model 
C is obtained from the extrapolation of the light distribution at 
latitudes larger than 5°. 

The mean durations predicted by Grenacher a/. 1999 
(model D) take into account the EROS efficiencies. Therefore, 
we can directly compare the predictions with our (uncorrected) 
values. This model predicts longer durations than observed, 
but here again it seems that the predictions could be adjusted, as 
they are sensitive to the value of the minimal lens mass. 

9. Guidelines for further interpretation 

As we always emphasized when presenting previous results to- 
wards the spiral arms, the fact that the distance distribution of the 
target sources is poorly known complicates the optical depth in- 
terpretation. Fig. 20 shows the expected optical depth as a func- 
tion of the Galactic longitude /, for diff'erent target distances, us- 
ing model 1 . This figure shows first that the impact of the bulge 
on the optical depth is significant only towards y Set. One also 
sees that the optical depth to distances smaller than 12 kpc is 
larger for / > 0, the near side of the bulge, because the number 
density of lenses is larger on this side; but on the contrary, the 
OGLE II collaboration (Udalski et al. 2000a) reported a larger 
rate for / < than for / > in their catalog of microlensing 
events in the Galactic bulge. This is due to the fact that the opti- 
cal depth averaged over all distances can be larger for / < 0, the 
far side of the bulge, because there are more distant sources on 
this side, with a larger optical depth. The observed asymmetry 
should then depend on the selection of the monitored sources. 
This illustrates the diff'erence between the optical depth up to 
a given distance, and the average optical depth over a set of 



stars located at various distances. Therefore, we provide here 
guidelines to interpret our optical depth values within a Galactic 
model framework. 

9.1. The concept of "catalog optical depth" 

The EROS measured optical depth is the average over the dis- 
tance distribution of the monitored sources. Establishing this dis- 
tance distribution through individual spectrophotometric mea- 
surements would require an enormous amount of complemen- 
tary observations. This leads us to define the concept of "catalog 
optical depth" Tcat^ which is relative to a particular catalog of 
monitored stars: Teat is defined as the fraction of stars of a given 
catalog that undergo a magnification A > 1.34. Our measured 
optical depth can be compared with the depth derived from a 
lens and source distribution model as follows: first, one has to 
generate a synthetic source catalog that matches our own cata- 
log, taking into account the EROS star detection efficiency (Fig. 
4); then one can use the generated source distance distribution to 
estimate the average optical depth and compare it with the mea- 
surements. This procedure is described in more detail below. 

9.2. Synthesizing a catalog that matches the EROS one 

Let Smodei(D, Ml, My - Mj) be the source number density pre- 
dicted by the model as a function of the distance D and of the 
absolute magnitude Mj and color My - Mj; let A/(D) and Ay(D) 
be the predicted absorptions in I and V. The pre-requisite for 
an optical depth interpretation is that ncatiL V - I), the den- 
sity per solid angle per apparent magnitude and color index 
of the synthesized source catalog, fits the observed one within 
the visible color-magnitude diagram. This density is related to 
Smodei(D, Mi, My - Mi) as follows: 

ricatil, ^ - ^) SmodeliD, Mi, My - Mi)€star(BEROS )D^dD (22) 

where estari^ERos) is our star detection eflficiency given in Fig. 4 
as a function of Beros = / + 0.6(y - /) apparent B magnitude in 
the EROS system. For given / and V -I apparent magnitude and 
color, Mi and My - Mi depend on the distance D as follows: 

Mi = I-2.5log(D/10kpc)-AiiD), (23) 
My -Mi = V - I - AyiD) -\- Ai(D). 

For the comparison with the observations, the color-magnitude 
diagrams nerosiL V - /) of our catalogs (Fig. 5) should be used. 
They provide the observed stellar density per unit solid angle 
in the eff'ectively monitored field, that is corrected for the dead 
regions of our CCDs, the overlap between fields and the blind 
regions around the brightest objects. 

We show here the color-magnitude diagram of a preliminary 
catalog obtained with this procedure (Fig. 21). This catalog has 
been produced for the qualitative understanding of the observed 
diagrams (Sect. 4.2, Fig. 5) and also of the color bias of the 
lensed population (Sect. 6.4.2, Fig. 11). For this simulation, we 
considered a stellar population identical to the one of the so- 
lar neighborhood (ESA 1997, Turon et al. 1995), distributed ac- 
cording to the disk and bulge mass distributions described in 
Sect. 8; the reddening was calculated assuming the following 
absorption law (Weingartner & Draine, 2001): 

A(V)/Nh = 5.3 X 10-^^cm^ A(/)/A^^ = 2.6 x IQ-^^cm^, (24) 

where Nh is the column-density of neutral hydrogen (atomic 
plus molecular), assumed to represent 4% of the mass density. 
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Fig. 21. Simulated color-magnitude diagram towards y Set (see text). 



The catalog was simulated using our star detection efficiency of 
Fig. 4. Fig. 1 1 (left) shows the source distance distribution of the 
catalog thus obtained. This catalog will be finalized in a forth- 
coming paper. 

9.3. Estimating the average optical depth 

Once the synthesis of a catalog is performed, one can compute 
the average optical depth over the synthesized catalog: 



^ model 



nmodel{D)Tmodel{D)D^ dD 



(25) 



^^n^odei{D)D^m 
where Tmodei(D) is the predicted optical depth at distance D and 

nmodel(D) = (26) 
JJ SmodeliD, Mi, My - M/) • estar{BEROs)^Mi d(My - M/) 

is the number density of stars (per solid angle and distance unit) 
located at distance D that enters the catalog. Beros = /+ 0.6(y- 
/) is obtained from Mj, My - Mj and D through equations (23). 

This procedure assumes that the average microlensing de- 
tection efficiency does not vary significantly with the distance of 
the sources. This is indeed the case because both the efficiency 
of our star detector and the average microlensing detection ef- 
ficiency depend on the apparent magnitudes. The microlensing 
detection efficiency averaged on distant detected stellar popula- 
tions is then similar to the one averaged on nearby detected pop- 
ulations because their apparent magnitude distributions are both 
peaked around the limiting magnitude. We have checked these 
features with simulations using our star detection efficiencies. 

The observed optical depths in Table 3 are relative to our 
entire catalogs. We recall here the average value for the 4 targets 
(12.9 X 10^ stars): 



^full catalog 



= 0.51+0.13x10" 



For easier comparisons, we extract from Fig. 16 the following 
optical depth 

f{Ic< 18.5) = 0.39 + .1 1 X 10-^ 

relative to the sub- set of 6.52 million stars brighter than = 
18.5, for which our catalog is close to being complete. 



10. Conclusions 

The microlensing event search of ER0S2 towards transparent 
windows of the spiral arms leads to optical depths that are con- 
sistent with a very simple Galactic model. The possibility of a 
long bar that was proposed in our paper II as a possible explana- 
tion of the observed optical depth and duration asymetries is not 
confirmed; indeed, the inclination of the bar was revised since 
the time of this publication, and the final optical depth measured 
towards y Set is also smaller (but statistically compatible). A 
more complete interpretation that would take into account the 
distance distribution of the monitored sources needs a model that 
allows one to synthesize the EROS catalogs of these sources. 
With such a model, one would also be able to make use of the 
event duration distributions. 

The VISTA project with its wide field infrared camera ap- 
pears to off'er an excellent opportunity to improve our knowl- 
edge of the microlensing towards the Galactic plane. The in- 
frared light will allow the observers to monitor stars through 
dust, making them free of the transparent windows that were 
Hmiting the EROS fields. 
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Fig. 22. Finding charts (84"x84", N-up, E-left) and magnification of the candidates as a function of the Heliocentric Julian Day (HJD-2, 450, 000). 
In the case of multi-detected events, we display only the best sampled light curve. 
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Fig. 23. Light curves and finding charts (continued) 
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Fig. 24. Light curves and finding charts (continued) 
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Fig. 25. Light curves and finding charts (continued) 



Mao S., Stefano R.D., 1995, ApJ, 440, 22 

MoUerach S., Roulet E., 2002, Gravitational lensing and microlensing. World 
Scientific 

OGLE web page http://www.astrouw.edu.pl/ ogle/ogle2/fields.html 
Paczynski B., 1986, ApJ, 304, 1 
Picaud S., Robin A.C., 2004, A&A, 428, 891 
Popowski, R, Griest, K., Thomas, C. L., et al, 2005, ApJ, 631, 879 
Rahal Y. R., 2003, Ph.D. thesis, Universite Paris 6 LAL-CNRS/IN2P3 report 
03-85 

Rahvar, S., 2004, MNRAS, 347, 213 

Sumi, T., Wozniak, R R., Udalski, A., et al. (OGLE coll.), 2006, ApJ, 636, 240 
Tisserand P., 2004, Ph.D. thesis, Universite de Nice-Sophya Antipolis CEA 

DAPNIA-04-09-T 
Tisserand R et a/. ,2007, A&A, 469, 387 
Thomas et al, 2005, ApJ, 631, 906 

Turon C, Requieme Y., Grenon M. et al, 1995, A&A 304, 82 
Udalski A., Szymahski M., Kaluzny J. et al. (OGLE Coll.), 1993, Act. Astr, 43, 
289 

Udalski A., Zebrun K., Szymahski M.,et al. 2000a, Acta. Astr 50, 1 
Udalski A., Szymahski M., Kubiak M.,et al. 2000b, Acta. Astr 50, 307 
Weingartner J.C., Draine B.T., 2001, ApJ, 548, 296 



